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SHOCK  HYDRODYNAMICS 
Mark  L.  Wilkins 

Lawrence  Radiation  Laboratory,  University  of  California 
Livermore,  California 

February  19,  1962 

ABSTRACT 

Thie  Hugoniot  relations,  which  express  conservation  of  mass,  momentum, 
and  energy  across  a  discontinuity,  are  derivedo  The  relations  are  applied 
to  a  perfect  gas  for  illustration,  and  graphical  methods  of  solving  shock 
problems  are  presented. 
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SHOCK  HYDRODYNAMICS 
Mark  L.  Wilkins 

Lawrence  Radiation  Laboratory,  University  of  California 
Livermore,  California 

February  19,  196Z 

INTRODUCTION 

This  article  is  divided  into  three  sections  that  are  not  necessarily 
dependent  upon  one  another.  The  first  is  the  derivation  of  the  Hugoniot 
relations  which  express  conservation  of  mass,  momentum,  and  energy  across 
a  discontinuity.  The  second  is  the  application  of  these  relations  to  a  perfect 
gas,  and  the  third  is  the  solution  of  shock  problems  by  graphical  methods. 

The  source  of  this  material  is  to  be  found  in  Hugoniot 's  original  paper 
in  the  Journal  de  L'Ecole  Poly  technique,  1889*  The  method  of  solving  shock 
problems  graphically  was  first  shown  to  me  by  Alan  Kaufman. 

In  this  analysis  a  shock  is  a  compression  wave  with  an  infinite  pressure 
gradient.  The  shock  width  is  small,  so  a  shock  wave  can  be  replaced  by  a 
surface  across  which  the  pressure,  density,  material  velocity,  and  internal 
energy  change  discontinuously. 

The  equations  of  state  used  are  in  the  form  of  pressure  as  a  function  of 
volume  and  internal  energy.  It  is  convenient  to  use  relative  volume  units,  or 
the  ratio  of  the  true  volume  to  its  original  volume.  An  equivalent  ratio  is  the 
ratio  of  the  original  density  to  the  present  density. 


1.  DERIVATION  OF  HUGONIOT  RELATIONS 

Given  a  fluid  in  an  initial  state  Eq,  p^,  Pg,  and  Ug  representing 
energy,  density,  pressure,  and  material  velocity,  respectively.  Consider  a 
shock  with  velocity  S  (with  respect  to  the  gas  velocity  in  front)  starting  from 
the  end  and  traveling  through  the  fluid,  changing  the  state  from  Eg,  pg,  Pg, 
Ug  to  E^,  p^,  P^,  U^.  In  a  time  £  a  length  L  =  St  will  have  been  swept  out 

(Fig.  la).  The  velocity  of  the  rear  surface  relative  to  the  fluid  is  (U^^  ^0^' 

therefore  the  rear  surface  will  have  been  displaced  (U^  -  Ug)t  during  the 
time  the  shock  has  traveled  the  length  L  =  St  (Fig.  lb). 
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1.  Conservation  of  Mass 

For  the  length  of  material  being  considered,  conservation  of  mass 
requires  the  mass  before  and  after  the  shock  has  passed  to  be  the  same.  The 
cross  section  is  considered  to  be  unity. 

StPo  =  [St  -  (Uj  -  Uj,)t]pj, 

^pQ  ^P^  "  I  ~ 

S{Po  -  Pi)  =  ”(Ui  -  Uq)p^, 


Z.  Conservation  of  Momentum 

Conservation  of  momentum  for  the  length  L  requires  that  the  net 
force  X  time  equal  the  change  in  momentum. 

(Pi  -  PqH  =  ppLUl  -  PqLUq 

Pi  -  Pq  " 

Substituting  (U^  -  U^)  from  equation  (I.l)  we  get 


Here  the  volumes  are  referred  to  the  density  p  ,  making  the  relative 
volume  Vq  equal  to  1.  Vq  is  carried  through  the  equations,  even  though  it  is 
1,  to  describe  the  general  case  where  the  volumes  are  referred  to  a  reference 
density  (p  -)  that  is  not  the  density  p^.  ahead  of  the  shock.  In  this  case  we 

would  have“V„  =  and  Vj  =  p^^p/pj. 
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Using  equations  (I.l)  and  (1.2)  to  eliminate  S  we  get  the  very  useful 
relation 

(I.2b) 


VoPo(«l  -  =  tPi  -  Po)(Vo  -  Vl)- 


3.  Conservation  of  Energy 

Conservation  of  energy  requires  that  the  net  work  on  the  mass  equal  the 
change  in  kinetic  and  internal  energy. 


(PfUi  -  PqUo)'  =  ^Po 


i  U, 


u 


E, 


where  and  Eq  represent  internal  energy  per  unit  mass. 

=  Spo[i(Ui  +  Uo><“i  ■  '^oQ  *  ■  ^oI 

(P.  -  P„)  _  _ 

'  B(Ui  +  Uq)(U^  -  Up)] 


Ui  -  Uo 


-  Uo 


+  VpoV 

U  1 1, 


from  equations  (I.l)  and  (I.Z). 

Ppi  -  +  “o’  +  v‘  -  V,-  <==i  -  '=o>Po'"o- 


"i-^o 


'Pi  +  Pq 


"i  -  Uo 


<Ul  -  Up)  =  v„  ■  V.  ■  Po>Po'^o’ 


p,  -1  p. 


(^1  ■'  Eo)Po'^0  =  '  2 . -  -  ^1>’ 

or,  including  PqVq  in  the  energy  units 

-1  -  -=0  =  ^  <Vo  -  Vl). 


(1-3) 


where  E  =  ^  units  of  the  reference  density). 

Equations  (I.l),  (1.2),  and  (1.3)  are  the  Hugoniot  relations  expressing 
conservation  of  mass,  momentum,  and  energy. 
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4.  Eliminating  Energy  Dependence 

For  a  given  equation  of  state  P  =  P(V,  E)  the  energy  dependence  may 
be  eliminated  by  the  third  Hugoniot  equation.  This  gives  the  Hugoniot  or 
pressure-volume  points  that  may  be  reached  by  a  shock  starting  from  a 
reference  state  (Fig.  2). 


P 


For  small  “  ^1  “  have  ~  ^2  ~  ^1  ~ 

=  P  dv,  or  the  adiabat  and  the  Hugoniot  coincide. 
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II.  APPLYING  THE  HUGONIOT  RELATIONS  TO  A  PERFECT  GAS 


1 .  Propagating  a  Uniform  Shock  into  a  Perfect  Gas  at  Rest 


Fig.  3. 


Equation  of  state: 


p  =  e 

V  MUL-16709 


Propagation  of  a  uniform  shock  into  a  perfect  gas  at  rest. 


Given  a  column  of  gas  at  rest,  at  one  end  a  piston  is  suddenly  given 
the  velocity  U  which  is  maintained  constant.  A  shock  S  travels  down  the 
column  changing  the  gas  from  the  state  subscript  zero  to  subscript  one.  We 
wish  to  find  the  new  state  variables  and  the  shock  speed  S. 


or 


1  ^0  PQ 

U  =  S{1  -  V^) 

(y  +  1)  -  V^(y 


.US 


p  =  p 
■^1  0 


-  1) 


V^(y  +  1)  -  (y  -  1) 


(Second  Hugoniot  equation),  (II.  1) 
(First  Hugoniot  equation),  (II. 2) 

(Third  Hugoniot  equation  (II-3) 

and  the  equation  of  state). 


2  +  (1  -  V^)(y  -  1) 
2  -  (1  -  V^)(y  +  1) 


Eliminating  S  from  equations  (II.  1)  and  (II. Z)  gives 


^1  "  ^0  +  TT 
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and  substituting  in  equation  (II. 3),  we  get 

,2 


(l  -  V,)-  . 


2  (Y+l)PoU‘ 


PQ 


2yPo  ^  ''"0 

Replacing  (1  -  V,)  by  U/S  from  equation  (II. 2)  we  get 


q2  y  +  1 

b  ^ 


US 


Pq 


=  0. 


(II.4) 


(II- 5) 


This  equation  gives  the  shock  speed  S  when  the  piston  velocity  and  the 
state  ahead  of  the  shock  are  known.  The  roots  of  the  equation  are  always 
real,  one  positive  and  the  other  negative  corresponding  to  whether  the  gas  is 
to  the  left  or  right  of  the  piston  that  is  compressing  it. 

The  positive  part  of  equation  (II.  5)  is 


_  Y  +  1 


U  -t- 


2 

+ 


PQ 


(II.6) 


The  speed  of  sound  ahead  of  the  shock  is  Cq  =  Pq-  Hence 


-  Y±J_ 


U  + 


fY  + 


+  c/. 


The  shock  speed  S  increases  with  U.  When  U  is  zero  S  -  Cq.  For 


large  U,  wliere  Cq/U  is  negligible,  S  - 


+  1 


U. 


P  -  P 

From  equation  (II.  1)  we  have  S  =  ---  •  Replacing  S  by  this  value 


Po 


in  equation  (II.  5)  we  get  an  expression  for  the  pressure  P^. 

(^1  -  Po)'  -  -  Po'  -  VPoPo"'  = 

Solving  equation  (11.7)  for  the  positive  root  gives 


(II. 7) 


Pi  =Po  +  H^  Oo'"' "  foP 

When  U  is  large  enough  that  Cq/U  is  negligible. 


+  c/. 


(II.8) 


^0  = 


=  (-4-0 


PoU". 
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The  relative  volume  behind  the  shock  is  obtained  by  eliminating  S 


from  equations  (II. 2)  and  (II. 6). 


Vi  =  l- 


u 


When  Cq/U  is  negligible, 


(H.9) 


V  =  T 1 
Y  +  1  ’ 


which  is  the  minimum  relative  volume  that  a  single  shock  can  produce. 


2.  Graphical  Representation 

To  get  a  clearer  view  of  these  relations  we  shall  represent  the 
Hugoniot  and  adiabat  graphically.  Taking  as  coordinates 
(remembering  that  Vj  is  the  relative  volume  behind  the  shock  where  the 
relative  volume  ahead  of  the  shock  is  1)  and  using  equations  (II. 3)  for  the 
Hugoniot  we  get  the  curves  shown  in  Fig.  4. 


Fig.  4.  Hugoniot  and  adiabat  starting  from  same  point. 
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The  Hugoniot  starting  from  point  A  has  the  asymptote  -  "I-  while  the 
adiabat  has  the  ordinate  as  an  asymptote.  By  differentiating  the  Hugoniot  and 
the  adiabat  two  times  with  respect  to  V we  get 


d(Pl/Po) 

r" 

d{Pi/Po) 

dVj 

A  (Hugoniot)  1 

_  J 

A  (adiabat) 


d^(Pj/PQ) 

d^P/Po 

J 

A  (Hugoniot) 

A  (adiabat) 


=  y(y  +  1). 


which  shows  that  the  Hugoniot  and  adiabat  have  at  point  A  the  same  tangent 
and  curvature,  as  was  pointed  out  in  section  I. 


3.  Reflection  of  a  Uniform  Shock 

In  the  preceding  analysis  we  have  considered  the  column  of  gas  to  be 
infinite  in  length.  Now,  we  shall  suppose  that  it  is  terminated  by  a  '^stone 
walH^  where  the  velocity  is  always  zero. 

When  the  shock  S  from  the  piston,  equation  (II. 6),  reaches  the  stone 
wall  a  reflection  is  produced.  That  is,  a  new  shock  is  formed  which 
travels  back  toward  the  piston  and  changes  the  gas  velocity  from  to  that 
of  the  stone  wall  or  zero.  This  shock  reaches  the  piston  and  a  new  shock 
S^,  analogous  to  S,  is  formed,  etc. 

We  wish  to  find  the  values  of  P  and  V  behind  a  shock  S,  where  for 

n  n 

odd  n  the  shock  is  traveling  from  the  piston  toward  the  stone  wall  and  for 
even  n  the  shock  is  traveling  from  the  stone  wall  toward  the  piston. 


S 

n 


(U  -  U  ,)V  , 

'  n-1  n-1 


n 


V  ,  -  V 
n- 1  n 


(II.IO) 


P 

n 


-  P 


n-1 


p  ,S  (U  -  U  J. 

^n-1  n'  n  n-P 


n'  n 


(II.  11) 


[Note  that  equations  (11.10)  and  (II.  11)  are  generalized  versions  of  equations 
(I.l)  and  (1.2),  respectively.]  The  quantity  (U^  -  j^)  alternates  from  (+) 

the  piston  velocity  to  (-)  the  piston  velocity  U.  For  a  shock  process, 

P^  >  P^  so  S  changes  sign  corresponding  to  its  direction.  Eliminating 
S  from  equations  (II.IO)  and  (II.ll)  and  taking  account  of  the  fact  that 

Pn-l  =  fo/Vr 
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V 

n 


=  V 


n-1 


P 

n 


PO 


u 


2 


-  P 


n-1 


(11.12) 


P 

n 


=  P 


n-1 


+ 


V  +  1 
4 


+ 


+  C 


n-1 


(11.13) 


[Equation  (11.13)  is  derived  from  equation  (II. 8).] 

The  solution  of  the  problem  is  complete  since  all  the  quantities  of  index 


n  can  be  gotten  from  the  quantities  of  index  n-1. 


4.  Application 

We  wish  to  apply  equations  (11.12)  and  (11.13)  to  get  the  conditions  behind 
the  first  reflected  shock  from  the  stone  wall.  Considering  the  initial  state  of 
the  gas  to  be  =  1,  P^  =  0,  and  density  p^,  we  get  from  equations  (11.12) 
and  (11.13)  the  state  behind  the  incoming  shock: 


=  V  ~ 

1  y  +  1 


+  1 


PQ 


P 


y  +  1 
Y  -  1 


pQ- 


Substituting  these  values  back  in  equations  (11.12)  and  (11.13),  we  get  the 
state  behind  the  shock  reflected  from  the  wall: 


III.  GRAPHICAL  METHODS  OF  SOLUTION 


Problem  1 

A  plate  Pj^  originally  traveling  to  the  right  with  velocity  Uq  strikes  a 
target  which  is  at  rest  (see  Fig.  5).  After  collision,  the  velocity  Uq  of  the 
plate  is  U^.  A  shock  S^  is  transmitted  to  the  right  into  the  target  and  a 
shock  S_,  is  transmitted  to  the  left  into  the  plate.  The  shock  S_  changes  the 

z  ^ 
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Fig,  5.  Plate  P,  traveling  to  the  right  at  velocity 
strikes  stationary  target,  creating  shocks  and  S  in 
plate  and  target  as  shown.  Shock  in  target  traveds  to 
interface  with  second  plate,  P^,  and  creates  transmitted 
and  reflected  shocks  S^-  and  S  ,  respectively. 
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target  material  from  state  P^,  E^,  to  state  P^,  E^,  P3  •  Similarly  the 

shock  into  the  plate  changes  the  plate  material  from  state  Pq,  Eq,  p^  to 

state  P  ,  E,,  p_,.  The  pressures  at  the  interface,  P  and  P  ,  will  be  equal. 
[See  equations  (III.l).] 


Equations  (III.l) 


For  the  plate  P^  For  the  target 


From  Hugoniot  mass  and  momentum  conservation  equations: 

Po(«3  -  =  <^2  -  -  ''z'l  Pl(”3  -  "1)^  "  <^3  -  -  ''3> 


Equation  of  state: 

^2  ^  ^plate^^2’  ^2^ 


^3  ■  ^target^^3’  ^3^ 


Third  Hugoniot  equation: 


^2  ^0 

^2  -  ^0  - 


0 


^2) 


E, 


^1  = 


P3  +  P 


^  (V,  -  V3) 


The  third  Hugoniot  equation  can  be  used  to  eliminate  E  from  the 
equation  of  state  to  give  the  "Hugoniot  equation"  for  the  material  where  P  is 
a  function  of  V  alone.  In  this  problem  we  will  consider  the  pressure  and 
energy  for  the  plate  and  target  to  be  zero  before  the  collision. 

Rewriting  the  equations  (III.l)  with  Hugoniots  for  the  equations  of  state 
we  have: 


Equations  (III. 2) 


For  the  target 


Pl\U3-Uj  =P2(Vi-V3) 

fv  1 

target 


P3  =  «.ar„t<''3' 


For  the  plate 


OoK  -  “0)"  '  P2l'^0  -  ''2' 


^2  =  “plat.(''2l 
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We  are  given  the  conditions  before  the  collision:  Uq,  Vq,  and  for 
the  plate,  and  V^,  and  for  the  target.  Since  there  are 

four  unknowns  —  P^,  V^,  V^,  and  —  and  we  have  four  equations.  The  four 
equations  are  not  readily  solved  algebraically,  so  a  graphical  method  will 
be  used.  Curves  of  P  vs  AU  are  made  for  each  material  by  assigning  values 
to  and  V^,  respectively,  in  the  equations  (III. 2)  and  calculating  the  resulting 
change  in  velocity  AU.  These  curves  (Fig.  6)  then  represent  all  of  the  possible 
pressure- velocity  states  each  material  could  have  as  a  result  of  a  strong 
shock  starting  from  an  initial  state  P  =  0  and  E  =  0. 

The  change  in  velocity,  AU,  may  be  positive  or  negative:  the  sign  is 
negative  if  the  material  is  slowed  down  from  a  given  velocity  by  hitting  an  ^ 
object  and  the  sign  is  positive  if  the  material  is  speeded  up  by  an  object 
hitting  it.  The  problem  of  a  graphical  solution  of  a  plate  hitting  a  target  is  to 
get  the  two  P-vs-AU  curves  into  the  same  coordinate  system. 

If  we  consider  the  coordinate  frame  of  the  target  where  in  this  case 
U ^  =  0  and  U^  is  an  increasing  velocity  to  the  right,  then  in  this  frame  the 
plate  has  P  =  0  when  U^  =  Uq  a-nd  the  pressure  will  increase  as  U^  decreases 
from  this  point.  Hence  the  solution  to  the  problem  is  obtained  by  super¬ 
imposing  the  curves  as  shown  in  Fig.  7. 

The  volumes  and  found  by  substituting  the  values  of  P  and 

U^  at  the  intersection  of  the  two  curves  into  equations  (III. 2). 

In  practice  the  Hugoniot  curves  for  various  materials  are  made  up  in 
advance  in  the  form  of  P-vs-AU  curves  where  the  initial  state  is  considered 
to  be  U  =  0,  P  =  0,  V  =  1,  and  E  =  0.  Then  the  left  side  of  the  P-vs-AU 
curve  is  obtained  by  ^Hlipping’*  the  curve  about  the  ordinate  as  was  done  above. 

From  Fig.  7  it  can  be  seen  that  if  the  target  and  plate  were  of  the  same 
material,  the  value  of  U^  at  the  interface  would  be  l/Z  the  incoming  velocity 
Uq.  If  the  target  material  were  infinitely  stiff  as  a  ^’stone  wall,  where 
U^  =  0  for  all  P,  the  pressure  at  collision  would  be  the  value  of  P  for  the 
incoming  velocity  Uq.  If  the  material  were  infinitely  soft  or  a  ^Void,  with 
P  =  0  for  all  U^,  then  pressure  at  collision  would  be  zero  and  U^  would 
remain  equal  to  Uq. 


Problem  2 

Consider  a  second  plate  to  be  at  rest  to  the  right  of  the  target  in  problem 
1  (see  Fig.  5).  The  shock  traveling  through  the  target  meets  the  second 
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plate,  P^,  and  transmits  a  shock  into  it,  changing  the  state  of  material 
from  subscript  6  to  subscript  5  as  shown.  The  reflected  signal  transmitted 
into  the  target  may  be  a  shock  or  a  rarefaction;  however,  the  signal  in  the 
plate  will  always  be  a  shock.  The  velocity  at  the  interface  is  U^,  and  again 
the  pressures  at  the  interface,  and  P^,  are  equal. 

First  we  shall  assume  to  be  a  shock.  Writing  the  Hugoniot  equations 
as  before,  we  have  equations  (III. 3) 

Equations  (III. 3) 


For  the  target 


For  plate  Z 


Pi  (U4  -  U3)"  =  (P4  -  PjXVj  -  V^)  p^(u^  -  =  (Pj  -  P^)(V^  -  V3) 


P4  =  «.arget<''4> 


P  =  H  ,  ^  ^(V^) 

5  plate  2  5 


We  are  given  the  conditions  just  as  the  shock  reaches  the  interface: 
U3,  V^,  and  for  the  target,  and  U^,  V^,  and  for  plate  2.  Since 


P4  =  P^,  again  there  are  four  unknowns  —  P^,  V^, 


V^,  and  U  .  —  and  four 
5  4 


equations.  The  equations  are  solved  graphically  as  before,  except  here  the 
Hugoniot  equation  for  the  target  is  not  the  same  since  it  is  the  locus  of 
pressure- volume  points  starting  for  E^  /  0  and  P^  7^  0.  It  is  obtained  in  the 
same  manner  by  substituting  the  third  Hugoniot  equation  into  the  equation  of 


state 
P  =  P 


Q.e., 


substitute  E. 


E, 


(V^  -  V4)  in  the  equation  of  state 


target^^4’  ^4^  eliminate  ^4  I  •  This  Hugoniot  compares  with 

the  Hugoniot  from  E  =  0,  P  =  0,  V  =  1  as  shown  schematically  in  Fig.  8. 

The  target  Hugoniot  starting  from  and  P^  and  the  plate  Hugoniot  are  used 
to  construct  P-vs-AU  curves  by  using  equations  (III. 3)  in  the  same  way  as 
was  done  before. 

The  graphical  solution  consists  of  getting  the  target  P-vs-AU  curve 
into  the  coordinate  frame  of  the  plate  as  shown  in  Fig.  9- 

The  volumes  and  are  found  by  substituting  the  values  of  P^  and 


U4  at  the  intersection  of  the  two  curves  in  the  equations  (III. 3). 
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For  most  solid  metals  the  energy  dependence  in  the  equations  of  state 
is  small  so  that  over  limited  ranges  the  Hugoniot  starting  from  V^, 
in  Fig.  8  turns  out  to  be  very  close  to  the  principal  Hugoniot.  Also,  over 
limited  ranges  the  adiabat  expanding  from  a  point  on  the  Hugoniot  will  lie 
close  to  the  Hugoniot.  These  facts  greatly  simplify  and  increase  the  range 
of  problems  that  can  be  solved  in  this  manner.  (See  Appendix  for  an  analysis 
on  a  typical  equation  of  state  for  a  metal.) 

Because  of  the  small  energy  dependence  it  will  not  be  necessary  in 
problem  2  to  construct  the  P-vs-AU  curve  starting  with  the  conditions 
V  =  V^,  P  =  P^;  E  =  E^,  and  U  =  instead  we  may  use  the  P-vs-AU  curve 
that  was  obtained  from  the  state  P  =  0,  E  =  0,  V  =  1,  and  V  =  0.  We  still 
have  the  problem  of  getting  into  the  right  coordinate  frame. 

We  want  to  match  the  velocities  so  that  the  target  is  moving  at  velocity 
in  the  coordinate  frame  of  the  plate.  This  is  accomplished  by  flipping  the 
target  P-vs-U  curve  about  as  shown  in  Fig.  10. 


Fig.  10.  Graphical  solution  for  P  and  U  at  interface,  using  the  principal 
Hugoniots  for  plate  and  target. 
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From  Fig.  10  we  see  that  the  target  pressure  increased  from  to 
indicating  a  shock.  If  plate  2  had  been  the  same  material  as  the  target  the 
P-vs-AU  curve  would  have  been  the  same  as  the  target  curve  before  flipping 
and  no  change  in  pressures  would  have  occurred  at  the  interface.  If  the  plate 
2  curve  had  intersected  the  target  curve  to  the  right  of  P^  the  pressure  would 
have  dropped  in  the  target,  indicating  a  rarefaction.  If  plate  2  had  been  a 
void  (P  =  0  for  all  U)  the  pressure  P^  would  have  gone  to  zero  and  the 
velocity  from  to  Thus  the  front  surface  velocity  is  equal  to  twice 

the  particle  velocity  for  a  shock  traversing  the  material  and  reaching  a  free 
surface . 

In  equation-of- state  measurements  the  free  surface  velocity  is  meas¬ 
ured  for  an  unknown  material.  From  this  the  particle  velocity  is  known,  and 
if  the  shock  transit  time  is  measured  the  pressure  can  be  obtained  by  the 
Hugoniot  relation  P  =  p^US.  The  other  Hugoniot  relations  give  the  volume 
and  energy.  A  series  of  P,  V,  E  points  can  be  obtained  in  this  way  by 
repeating  the  experiment  with  various  input  velocities.  These  points  do  not 
allow  one  to  immediately  write  the  pressure,  volume,  energy  relations  or 
the  surface  that  represents  all  the  states  of  the  material  since  we  only  have 
the  experimental  points  of  a  line  on  this  surface.  To  write  the  equation  of 
state,  assumptions  are  made  concerning  the  energy  dependence  and  a  P,  V, 

E  relation  is  written  satisfying  the  experimental  Hugoniot  points. 

Problem  3 

Consider  a  length  of  high  explosive  (H^E.)  that  has  been  detonated,  where 
the  detonation  front  has  just  reached  the  surface  of  an  adjoining  plate  at  rest. 

A  shock  is  transmitted  into  the  plate.  The  problem  is  to  find  the  pressure^ 
density,  and  internal  energy  behind  the  shock.  (See  Fig.  11.) 

The  three  conservation  laws  must  still  hold.  The  problem  is  similar 
to  problem  2,  but  here  the  detonation  front  replaces  the  shock  The 

Hugoniot  curve  for  the  He  E  is  constructed  starting  from  the  conditions  at  the 
detonation  front,  called  the  Chapman- Jouguet  point  (C-J),  and  the  H.E.  equa¬ 
tion  of  state,  all  considered  to  be  known  in  advance.  The  procedure  is 
identical  to  problem  2,  with  the  subscript  3  being  replaced  by  the  values  at 
the  C-J  point.  The  P-vs-U  curves  for  both  materials  are  superimposed  in 
Fig.  12  (just  as  was  done  in  Fig.  9)j  and  the  intersection  gives  the  values  of 
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(A. 2a) 


The  underlined  terms  in  equation  (A. 2a)  are  the  same  as  the  principal 
Hugoniot.  Since  ^  and  c  is  usually  between  1  and  2,  equation  (A. 2a) 
is  very  nearly  equal  to  the  principal  Hugoniot  except  for  high  order  terms. 


APPENDIX  B 

To  show  that  the  front  surface  velocity  equals  twice  the  particle  velocity 
(Ufs  =  2U  )  for  a  material  with  a  Hooke's  law  equations  of  state, 

P  =  a(p/p^Q  -  1). 

=  —  =  — 

9p  Pq 

where 

a  =  constant, 

Pq  =  reference  density, 

c  =  sound  speed. 

The  Riemann  invariant,  cr,  for  the  hydrodynamic  equations  of  motion  is 

O’  =  ^  c-^.  (See  R.  Courant  and  K.  O.  Fredricks,  Supersonic  Flow  and  Shock 

Waves,  Inter  science,  New  York,  1948.) 

Since  c  is  a  constant  for  this  equation  of  state  we  have 

or  =  c  \  =  c  In  =  c  In  n, 

J  P  pQ 

^  =  p/  Pq- 

For  a  shock  traveling  through  undisturbed  material  we  have 


where 
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In  the  above  derivation  the  equation  of  state  v/as  assumed  to  be  linear 
in  T|.  A  more  critical  analysis,  using  a  nonlinear  equation  of  state,  v/ill  show 
the  preceding  result  is  valid  over  a  large  range  of  pressure  and  compression. 
Consequently  the  assumption  of  using  the  principal  Hugoniot  in  solving  shock 
interaction  problems  is  reasonable. 
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This  report  was  prepared  as  an  account  of  Government 
sponsored  work.  Neither  the  United  States,  nor  the  Com¬ 
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